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We discuss a phase diagram of two-dimentional U{1) x f/(l) superconductor in the field theoretic 
formalizm of Ref. [17]. In particular we discuss that when penetration length is short the system 
exhibit a quasi-neutral quasi-superfluid state which is a state when quasi-long range order sets in 
only in phase difference while individually the phases are disordered. 



INTRODUCTION 



TWO-GAP GINZBURG-LANDAU MODEL AND 
ITS REPRESENTATION IN GAUGE-INVARIANT 
VARIABLES 



Berezinskii-Kosterlitz-Thouless (BKT) transitions and 
phase diagrams of planar superfluids and superconduc- 
tors is a subject of wide interest and intensive studies [1]- 
[10]. The present study is motivated by condensed matter 
systems with multiple coexistent condensates [11]- [14]. 
In two-band superconductors a Fermi surface consists 
of several disconnected parts which gives rise to several 
types of carriers (which are "living" on different parts of 
the Fermi surface) and correspondingly to multiple gaps 
[11]-[13]. In known two-band superconductors the two 
types of carriers are not independently conserved and 
the U{1) X U{1) symmetry is softly or strongly broken to 
U{1). The two-gap model with the exact U{\) x U{1) 
symmetry has been discussed in the theoretical stud- 
ies of liquid metallic hydrogen [14]. In such a system 
(which appears being close to realization in high pres- 
sure experiments) there should be superconductivity of 
electron-electron and proton-proton Cooper pairs which 
gives rise to the U{\) x U{1) symmetry. Besides that the 
two-gap superconductivity is indeed a question of an ab- 
stract academic interest due to this system albeit being 
formally very simple has deep physics and exhibit numer- 
ous interesting counter-intuitive phenomena (see also re- 
mark [15]). We also mention that certain formal aspects 
of a BKT transition in a square-lattice planar two-gap 
Abelian Higgs model were earlier also a subject of study 
in context of the lattice theory which was considered in 
connection with spin-charge separated superconductivity 
[10]. However, physically, the phase diagram of two-gap 
system is principally different from the spin-charge sep- 
arated superconductor. In the papers by Rodriguez [10] 
it was analyzed a two-component model which was sug- 
gested to be relevant for spin-charge separated supercon- 
ductor in a presence of a coupling to statistical gauge 
field. The main principal distinction of the system [10] 
and a two-gap superconductor is that in a spin-charge 
separated superconductor there exist two vector poten- 
tials which correspond to statistical and ordinary gauge 
fields [16], [10]. 



We begin with a brief outlining of basic properties of 
the Ginzburg-Landau functional for a two-gap (distin- 
guished by index a — 1,2) planar superconductor: 
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where ^'r 



and V{\^^\^) 



^l^'cj and rj is a characteristic of interband Josephson 
coupling strength [12]. In a general case (1) may include 
other potential terms and mixed gradients terms which 
we drop in order have discussion in the simplest form. 
Similar models are also discussed in particle physics [19]. 
The situation is quite complex in a charged system (1) 
since the condensates are coupled by the field A. How- 
ever the variables in (1) can actually be separated in the 
London limit in a simply-connected space. In [13] it was 
shown that the model (1) is exactly equivalent to an ex- 
tended version of Faddeev 0(3) nonlinear cr-model [20] 
(which is also relevant in high energy physics [21] and for 
triplet superconductors [22]). The version of this model 
discussed in [13] consists of a three-component unit vec- 
tor n in interaction with a vector field C and a density- 
related variable p: 
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and 



V = A + Bns + Cnl (3) 
A — p^[4cim^ -|- Ac2m\ — hiuii — &2'7i2] 
B — p^[8c2m2 — Sciml — &2'7^2 + bimi] 
C = Ap^[ciml + C2ml]. 
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A position of the unit vector n on the sphere .5^ can be 
characterized by two angles as foUows: 

n = (sin(6') cos (7„), sin (6*) sin(7„), cos(^)). (4) 

The variables of (1) and (2) are related in the following 
way [13]: 
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p^Km = ?7[*i*2 + *2*i]> where K = 2r)y/mTm^ is 



the Josephson term. The field C is directly related to 
supercurrent 
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The mass of the field C is the manifestation of the Meiss- 
ner effect [13], the corresponding magnetic field penetra- 
tion length is: 
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The potential term V in (2) breaks exact 0(3) symmetry 
associated with n to 0(2) and defines the ground state 
for ns (we denote it by ns), which corresponds to uniform 
density of condensates [13] : 
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where A^i,2 stands the for average concentrations of 
Cooper pairs < |*i,2|^ >• The term 

p^Kni breaks the remaining 0(2) symmetry. Thus 
the ground state of n is a point where cos7„ = — 1 if 
r? > or cos7„ = 1 if < on a circle specified by 
the condition = n^. In case when 77 w the system 
has an approximate U{1) symmetry associated with 7„. 
With this we conclude the brief reminder of the basic 
properties of the model (1), more details can be found in 
[13, 17]. 

Indeed we can study the phase diagram of a two-gap 
superconductor cither using the variables ni and C or the 
initial variables |\E'a|e~*"^°. The former representation is 
more compact and convenient, while the later variables 
are more traditional so below we discuss the phase dia- 
gram of the model in terms of both representations. Ad- 
ditional motivation for it is that besides superconductiv- 
ity the model (2) is of independent interest being relevant 
in particle physics [21]. 



VORTICES IN TWO-GAP SUPERCONDUCTOR 

In order to understand BKT transitions in a two-gap 
system one should understand nature of topological ex- 
citations. The two-gap system allows vortices with frac- 
tional flux [17] (see also the remark [18]). In this section 
we discuss the physical origin of the vortices in the models 
(1) and (2) (somewhat extending the discussion in [17]). 
Before we discuss vortices with fractional flux, we would 
like to rcnnind briefly the origin of flux quantization in 
ordinary superconductors. First let us however consider 
a topological defect in a neutral U{1) system like *He 
which is described by a complex scalar field |^'|e"^ with 
free energy 
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Such a system enjoys topological excitations of the 

map (when the phase </> changes 27rn around the core), 
in the form of vortices which have logarithmically diver- 
gent energy per unit length [23] . 

Let us now consider a superconductor. That is a 
charged system, so we have a coupling to the vector po- 
tential A 

Let us consider a vortex where the phase changes 27r 
around the core (27m winding of the phase is a topo- 
logical requirement which needed for singlevaluedness of 
the order parameter). Then in contrast to neutral case, 
away from the vortex core we can minimize energy by 
compensating gradients of phase in the first term of (11) 
by a corresponding configuration of vector potential A, 
which, as shown by Abrikosov [24], makes this defect be- 
ing of finite energy per unit length. The corresponding 
solution for the vector potential for a vortex character- 
ized by 27r phase change around its core is following: 



A = 



r X e. 



|A(r)| 



\A{r)\ = --J—KJeJ—r 
er V m \ V m , 
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where r measures distance from the core and is a unit 
vector pointing along the core. 

As it was discussed in [24] the above solution for A 
corresponds to the situation when a vortex carries a mag- 
netic fiux quantum $0 = 27r/e. Thus the flux quantiza- 
tion in the ordinary superconductor is an effect rooted in 
topological and energetical considerations. 

Let us examine now the two-gap model (1) in case of 
1] = and show that the two-gap Abelian Higgs model 
can not generate magnetic flux which would compensate 
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gradients of both order parameters in case of a general 
topological defect. 

Let us consider a vortex in (1) such that the phase 0i 
has 27r winding around core while (f>2 is constant. Then if 
a system is "trying" to compensate gradient of (pi in the 
first term by a corresponding configuration of A (that is, 
by generating flux) at the same time the flcld A is in- 
ducing divergence in the second term. So, by generating 
exactly one flux quantum the system may make the first 
term being finite but it would induce a divergence in the 
second term. On the other hand, by decreasing flux, the 
system may minimize divergence in the second term but 
then the gradient in the first term is not compensated 
by A identically and also becomes divergent. Thus there 
are no finite energy solutions for this type of vortices. 
However, still there is a solution with a minimal energy 
[17] (a reminder of derivation is given below), which in- 
deed should be weakly divergent like in a neutral U{1) 
system. Apparently also in the case where the phase ipi 
has 27r winding around core while 4>2 has a — 27r winding, 
the system may minimize one of the gradient terms by 
generating flux but at the same time it would increase 
divergence in the second term. That is because to mini- 
mize the configuration with winding 2tt one should have 
magnetic flux going in one direction while for the case 
— 27r in opposite direction. In the case where the phase 
(f)i has 27rn winding around core and (j)2 also has 2nn 
winding the system can generate flux which would com- 
pensate both gradients [17]. Only this type of topological 
defects in this system has finite energy. 

Let us now describe formally the vortices with 
fractional flux following [17]. In the London limit 
(|^'a|=const) the eqs. (1), (2) become: 
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Eq. (14) can also be expressed in traditional variables: 
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which is a rearrangement of the variables in (1) with no 
apprioximations involved. The first term, in the case 
when 77 = 0, describes a neutral boson (a coupling to 
A is eliminated when we extract the difference of gauge 



invariant phase gradients) associated with n (with 0(3) 
symmetry of n, in the London limit being broken down 
to 0(2) by mass terms for n^). The mass term for rii 
breaks the remaining 0(2) symmetry. The second term 
is the same as the second term in (14), together with 
third term 6^/2 = (32e2)-i[ViQ - VjCjj^ it describes 
the massive vector field C. 

Let us observe that the models (1) and (2) have four 
characteristic length scales - coherence lengths of conden- 
sates, magnetic field penetration length A and also there 
is a length scale associated with intrinsic Joscphson effect 
(which is inverse mass for m). Below we discuss a sys- 
tem where both coherence lengths are of the same order 
of magnitude and where the characteristic length scale 
associated with inverse mass for m is much larger than 
magnetic field penetration length A in type-II limit and 
much larger than coherence lengths in typc-I limit (which 
amounts to a negligibly small influence of the Josephson 
coupling at length scales much smaller than the inverse 
mass for rii). The discussion of the phase diagram of the 
system in the regime of large r] will be presented in an- 
other publication [25]. We stress that in a rigorous sence 
we should have 77 = in order to speak about true BKT 
transitions, however when r] is nonzero but sufficiently 
small we can speak about some finite-size crossovers. We 
also note that in principle it is possible to have physi- 
cal systems where 77 is exactly zero - such situation ap- 
pears e.g. in theoretical studies of liquid metallic hydro- 
gen (which appears being close to being realized in high 
pressure experiments) which should allow superconduc- 
tivity of protonic and electronic Cooper pairs [14]. Since 
electrons can not become protons the condensates are 
independently conserved and we have exact U{1) x U{1) 
symmetry. So in the zevo-T] case, which we consider in 
this paper, the system (1), (2) has a neutral U{1) sym- 
2mctry associated with the variable 7„ = — (t)2- 

Vortices in this system, characterized by the following 
, .phase change around a core A0i = 2iTki; A02 = 2iTk2, in 
^ 'what follows we denote for brevity by (fci, k2). From (14) 
and (15) one can observe that for vortices characterized 
byA(^!>i-h?!>2) = f^d/[V((/)i-K(/)2)] =47rm; A((/)i-(/)2) = 
0, (where we integrate over a closed curve a around a 
vortex core) the first term in (14) and (15) is identically 
zero and such vortices is a analogue of m-flux quanta 
Abrikosov vortices in an ordinary superconductor charac- 
terized by = I 

m \ 

that if both phases 
then a vortex carries one quantum of magnetic flux. 

Let us now consider the case when A^i = 27r and 

A(/)2 = [a vortex (1,0)]. First of all such a vor- 
tex features a neutral superflow characterized by a 27r 
gain in the variable 7„. As follows from (14), (15), 
the vortex in this respect is equivalent to a vortex in 
a neutral superfluid with superfluid stiffness \p^ s\v? 9 = 



I* 



l*2l' 
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+ ^ . Besides that, as follows from (14) 

and (15), a topological defect (1, 0) is necessarily accom- 
panied by a nontrivial configuration of the charged field 
C. That is, for such a topological defect the second 



term in (14) becomes 
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Thus, 



•}i — eA 

this term is nonvanishing for such a vortex configura- 
tion, which means that this vortex besides neutral vor- 
ticity also carries magnetic field. Here we stress that a 
minimal energy solution for a vortex with a given wind- 
ing number is a solution where the second term in (14) is 
made finite by a corresponding configuration of A. Lets 
us calculate the magnetic flux carried by such a vortex. 
The supercurrent around the core of this vortex is: 
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Lets us now integrate this expression over a closed curve 
a situated at a distance larger than A from the vortex 
core. Indeed at a distance much larger than the pen- 
etration length the supercurrent J, or equivalently the 
massive field C, vanishes. Thus we arrive at the follow- 
ing equation: 
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where <^ = §^ A.d\ is the magnetic flux carried by our vor- 
tex and $0 = 27r/e is the standard magnetic flux quan- 
tum. Obtaining equations of motions for C from (14) we 
can find a solution for asymptotic behavior of the field 
C for this vortex. At a distance larger than coherence 
length ^ = max[^\, ^2] from the core the field C behaves 
as: 



$0 



(18) 



where r is the distance from the core. So this vortex has 
the following structure: 

• r < ^ vortex core; 

• ^ < r < A this region features neutral superflow 
associated with the gradients of the variable 7„ as 
well as charged supercurrent associated with the 
field C. 

• r > A this region features neutral superflow (which 
like in a vortex in any neutral system vanishes as 
1/r away from the core), whereas at the distance 
r > A the magnetic field and the field C vanish 
exponentially according to (18). 

We can also derive the energy per unit length of this 
vortex which is: 



E = E^ + 
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where Ec is the core energy and R is the sample size. The 
second term in (19) is the energy of charged current and 
magnetic field which comes from the second and third 
terms in (15). The third term in (19) is the kinteic energy 
of neutral current which comes from the first term in (15). 
In an infinite sample, the energy per unit length (19), of 
this vortex in logarithmically divergent (like in a neutral 
superfluid) as a consequence of the existence of vorticity 
in the field of the massless neutral boson. When r] ^ 
vortex energy diverges linearly. 



We would like to stress that since the two U{1) symme- 
tries are not independent in this system but are coupled 
by vector potential the system actually behaves like a sys- 
tem with a,U{l) neutral and a J7(l) gauged symmetries 
as it is seen from separation of variables (14), (15). The 
physical origin of the massless neutral U{1) boson and 
flux fractionalisation in the two-gap system can be un- 
derstood also from the following considerations: We can 
write expressions for individual supercurrents in the Lon- 
don limit: j„ - ;^|*apV(/)„ - ^\M>^\^A, (a = 1,2) 
from where it explicitly follows that even if there are gra- 
dients of while there are no gradients of <j)2 nonetheless 
coupling by A induces current ]2 in the condensate '^2- 
For the vortex (1,0) such current partially compensates 
magnetic flux induced hy'^i and this compensation leads 
to existence of fractional flux and an effective neutral su- 
perflow in the system. The fractional flux follows from 
the condition that at a distance larger tham A from the 
vortex core the resulting charged current J = ji -I- j2 
should vanish. The physical meaning of the neutral su- 
perflow is the circulation of two sorts of charged Cooper 
pairs in opposite directions where the charged individual 
currents Jq compensate magnetic field induced by each 
other. 



We also remark that a straightforward inspection of (2) 
shows that this model could allow neutral vortices asso- 
ciated with the neutral 0{2) boson without a nontrivial 
configuration of the field C for any values of fi^- However 
such solutions (e.g. vortices in the case 713 7^ character- 
ized by A7„ = 27r; C = 0) are unphysical because these 
vortices do not satisfy the condition that changes by 
2'Kki around the vortex core (as follows from (14)). Thus, 
while the original model (1) and the extended Faddeev 
model (2) have the same number of degrees of freedom so 
that the fields n and C arc dynamically independent by 
construction, the mapping incurs a constraint on topo- 
logical defects in n and C: when we go around the vortex 
core the two phases (t>i should have 2'Kni gains (with rij 
being necessary integer). 
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BEREZINSKII-KOSTERLITZ-THOULESS 
TRANSITIONS IN TWO-GAP SYSTEM 

Wc begin with a remark that as it is well known, for- 
mally, the Abelian Higgs model H = 2^|(iV-eA)*)p + 
+ 1^4 jjQgg jjQ^ exhibit the BKT transition due 
to interaction between Abrikosov vortices is of short 
range, being screened over the magnetic field penetra- 
tion length A [3] (see also a remark [26] and a review 
[1]). However when the penetration length is signifi- 
cantly large the system undergoes a "would be" BKT 
transition or crossover (see a review [1]) approximately 

at the temperature T"^^ = | — — '-. The equation 
for T^^'^ should be solved self-consistently with equa- 
tion for the temperature-dependent gap modulus |\1/(T)|. 
The gap modulus 1^1/ (T)| opens at a characteristic tem- 
perature T* > T^^^ , and albeit the phase is random at 
T* >T> T^^'^, however |*(T)| has roughly the mean- 
ing of the measure of the density of preformed Cooper 
pairs in the system [27]. The "would be" BKT tran- 
sition in a charged system is experimentally observable 
in extreme type-II superconductors [1, 28). This transi- 
tion/crossover is associated with onset of superconduc- 
tivity in an extreme type-II planar system. In contrast a 
system with a short penetration length does not exhibit 
planar superconductivity at nonzero temperature. 

Let us now discuss the phase diagram of planar U{l)x 
U{1) superconductor in the type-I limit where both co- 
herence lengths are larger or equal to magnetic field pen- 
etration length; 



Phase diagram of a planEir two-gap system in the 
type-I limit 

At finite temperature a planar system generates vor- 
tices and antivortices. The BKT transition manifests it- 
self as formation of bound vortex-antivortex pairs. In the 
regime A < ^ we can neglect the presence of the second 
and third terms in (15) while describing interaction of 
vortices. The second and third terms describe charged 
supcrcurrent around the vortex core, it can be seen from 
(14) where second and third terms correspond to the sec- 
ond and third terms in (15) and thus these terms describe 
a massive vector field C. The inverse mass of this field 
is A thus a solution of equations of motion for C should 
be localized at the length scale A and in the limit of 
small A can be neglected (sec also a discussion after eq. 
(18)). So, when the penetration length is very short then 
we can neglect interaction of vortices mediated by the 
charged current which is exponentially suppressed. On 
the other hand the neutral mode [described by the first 
term in (14), (15)] is not affected by smallness of penetra- 
tion length. Thus effectively our system is described by 
the gradient term of the composite neutral mode which 



is (we remind that we consider the case of the negligible 
small rj): 
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We observe that the system has the following vortices 
with identical configuration of the neutral superfiow (1, 0) 
and (0, —1) and antivortices (—1, 0) and (0, 1). A neutral 
superfiow gives rise to a long-range logarithmic interac- 
tion between vortices and antivortices. So, these vortices 
undergo a genuine BKT transition at the temperature 
which can be derived immediately from the first term in 
(14) and (15): 



rBKT 



p2 sin^ e = 



2 _|_ |1'2P(Tf/^) 
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which, in the case of two-gap superconductor, should be 

solved sclf-consistently with the equations for the gap 
modules |\E'i,2(T')| [27]. We stress that in a rigorous sence 
we should have 77 = in order to speak about BKT tran- 
sitions, when 7] is nonzero then fractional vortices of the 
type (1,0) interact linearly at length scales larger than 
inverse mass for m and thus do not exhibit true BKT 
transition. However when r] is sufficiently small wc can 
speak about some finite-size effects which, under certain 
conditions may even be observable in an experiment. 

We should first emphasis that the transition (21) 
in type-I limit is associated with formation of vortex- 
antivortex pairs of the following types [(1,0) -|- (-1,0)], 
[(1,0) + (0,1)], [(-1,0) + (0,-1)], [(0,-1) + (0,1)] (due to 
e.g. the vortices (1,0) and (0,-1) have equivalent con- 
figuration of composite neutral superfiow in the limit 
A ^, and thus in Coulomb gas mapping represent iden- 
tical charge). We stress that the interaction of vortices 
with windings of phases belonging to different conden- 
sates like e.g. pairs [(1,0) + (0,1)] is mediated by neutral 
superfiow which has a composite nature. That is, this su- 
perfiow consists of particle current of both condensates 
(see a remark in conclusion of the previous section) . Be- 
sides that the system allows additional relevant compos- 
ite topological excitation: the one-flux-quantum vortices 
(aslo described in [17]) (1,1) and (-1,-1). These com- 
posite vortices are topologically equivalent to the pairs 
[(1,0) + (0,1)], [(-1,0) + (0,-1)]. The composite vor- 
tices (1,1) and (-1,-1) play an important role in this sys- 
tem: these vortices do not have neutral superfiow (first 
term in (15) is identically zero for such a vortex con- 
figuration) thus the interaction of these vortices is ex- 
ponentially screened, they have finite energy, and remain 
liberated at any nonzero temperature. Thus, due to exis- 
tence of the mentioned above four equivalent vortex pair- 
ing channels and the relevant composite vortices (1,1) 
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and (-1,-1), the phase transition (21) marks onset of the 
quasi-long-range order only in the variable 7n = — 4'2 
(described by first term in (14), (15)), while the variable 
7c = sin^(|)</)i+cos^(|)(/)2 associated with charged boson 
in (14), (15) remains disordered. In the system considered 
by Rodriguez [10] the counterpart of this transition is a 
transition in the physical electronic phase (the supercon- 
ducting transition) thus, the initially massless composite 
mode in [10] is actually coupled to gauge field like in an 
ordinary one gap superconductor and the phase transi- 
tion in the composite field is transformed in a "would 
be" BKT superconductive transition/crossover like in a 
one-gap Abclian Higgs model. In contrast, we stress the 
two-gap superconductor is a physical system which pos- 
sesses a genuine and potentially experimentally observ- 
able massless neutral boson associated with the compos- 
ite phase field 7„. Let us once more stress the meaning 
of the field 7„: 

• In a two-gap system variables can be separated into 
a phase diff'erence and a phase sum. Only the phase 
sum is coupled to vector potential. 

• The gradients of phase sum arc gauged away at the 
length scale A. (Just like in ordinary Abrikosov 
vortex a phase gradient is compensated by a vector 
potential at the penetration length scale). 

• When A is small the long- range interaction is me- 
diated only by gradients of phase difference which 
is not coupled to vector potential and thus is of 
long range and determines the energy of e.g. an 
(1,0)-|-(0,1) pair. When such pairing occurs the in- 
dividual phases are indeed disordered [in contrast 
to the case if pairing would be only of type (l,0)-|-(- 
1,0) and (0,l)-|-(0,-l)] and quasi-long-range order 
sets in only in phase difference. Away from a pair 
(1,0)-|-(0,1) the phase difference is not character- 
ized by a nontrivial winding number and its gradi- 
ents do not contribute to kinetic terms of the GL 
functional. At the same time the phase sum has a 
nontrivial winding but it is compensated by vector 
potential thus it gives only a finite contribution to 
the GL free energy which is not important for the 
BKT transition (similar like a finite energy contri- 
bution from vortex cores). So the system can be 
mapped onto a Coulomb gas in a standard way. 

We should stress that in two-gap superconductor (1) 

the transition in the field 7„ has little to do with super- 
conductivity since it is a transition in the composite neu- 
tral field. It marks onset of what should be rather called 
superfluidity in the field of the composite neutral boson. 
Also the important fact is that an elementary topological 
defect with neutral vorticity in the models (14), (15) also 
carries a fraction of magnetic flux quantum. Thus al- 
beit we have the phase ordering transition in the neutral 



field however the fractional flux which accompanies the 
vortices [17] allows a detection of this superfluid phase 
transition in a two-gap superconductor by standard ex- 
perimental technique used in studies of charged systems. 
That is, the vortex-antivortex interaction is mediated 
only by composite neutral superflow while the magnetic 
flux (at length scale A from the core) allows observa- 
tion of these defects in flux-noise measurements or mea- 
surements in applied dynamic magnetic field. Besides 
that because of the existence of neutral mode, the tran- 
sition (21) in e.g. fiux-noise measurements should man- 
ifest itself as a sharp well-defined transition which can 
not be observed in ordinary charged systems where one 
can see only a "would be" BKT transition, "washed out" 
by Meissner effect. So, a two-gap superconductor dis- 
plays unique potentially experimentally observable novel 
phase: in a type-I limit two-gap superconductor should 
exhibit quasi-superfluidity at the same time it should not 
exhibit true superconductivity at any nonzero tempera- 
ture. In the type-I limit the true superconductivity sets 
in at T — > because at zero temperature the dimensional 
reduction docs not work and the system is effectively 3- 
dimensional. In the intermediate case when penetration 
length is not negligibly small there could occur a washed 
out crossover to a superconducting phase at a nonzero 
characteristic temperature (T^/ >> Tc > 0). 

Let us now discuss the phase diagram of a two-gap 
system in the case of a large penetration length. 



Phcise diagram of a planar type-II two-gap system. 

When penetration length is much larger than coher- 
ence lengths, the system, again, possesses the following 
vortex excitations: (1,0), (-1,0), (0,1), (0,-1) which still 
all interact with each other by means of the composite 
neutral superflow. Also there exist composite vortices 
without neutral superflow (1,1) and (-1,-1). In contrast 
to the type-I limit, at larger A the system has two en- 
ergetically preferred vortex-antivortex pairing channels 
(instead of four equivalent vortex pairing channels) asso- 
ciated with binding of the following pairs [(1,0) -|- (-1,0)] 
and [(0,-1) + (0,1)]. It is because, as it follows from (14) 
and (15) the vortices (1,0) and (-1,0) have attractive in- 
teraction both due to neutral superflow and charged cur- 
rent, in contrast the vortices (1,0) and (0,-1) interact only 
due to their neutral modes [which is a composite mode, 
whereas the charged supercurrent which accompanies the 
vortices (1,0) and (0,-1) consists only of Cooper pairs of 
sort "1" in the case (1,0) and sort "2" in the case (0,-1) 
(see the concluding remark in section III)]. 

Let us illustrate the discussion with the vortex (1,0). 
The free energy density for this configuration becomes: 
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Prom this expression it is seen that the interaction of 

vortices (1,0) and (-1,0) at the length scales smaller than 
A is characterized by the effective stiffness 



J^ff = J" + , where J"(T) 
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(corrections due\ 
to finite A / 
(24) 

The contribution to the interaction between these vor- 
tices which is mediated by charged current (character- 
ized by the stiffness J°) should receive corrections due 
to charged current vanishes over the magnetic field pen- 
etration length A. Similarly, the effective stiffness for 
interaction between vortices (0,1) and (0,-1) is given by 
an analogous expression 



where 
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(25) 



(26) 



Thus the system in the regime A — > oo will have two 
transitions at the following temperatures: 



BKT\ 



-iBKT 



(27) 



This is a consequence of the fact that making A large, 
reduces four vortex pairing channels (which are equiv- 
alent in the type-I limit) to two non-equivalent ener- 
getically preferred pairing channels. We should em- 
phasis that the composite co-centered excitations (1,1) 
and (-1,-1) which do not carry superflow [and which are 
partially responsible for disordering the variable 7c = 
sin^(|)0i -I- cos^(|)^2 in the short penetration length 
limit] are irrelevant if the penetration length is large since 
the effective stiffness of interaction between vortices (1,1) 
and (-1,-1) is J = [|*ip/mi + |*2p/m2] (as it follows 
from (15); (14)) and thus when we can not assume that 
interaction is completely screened due to Meissner effect, 
then the large effective stiffness which characterizes in- 
teraction of these defects prevents existence of liberated 
vortices of this type at the temperatures (27). We also 
remark that in the type-II limit the interaction of vortex 
pairs of the type (1, 0) + (0, 1) is getting depleted at the 
length scales r < A as shown in appendix A. 



CONCLUSION 

We presented a discussion of the phase diagrams of 
planar U{1) x U{1) superconductors and predicted what 
could be a novel quasi-superfluid state realized in two-gap 
superconductor in a type-I limit. The physical meaning 
of the quasi-superfluid state in a two-component charged 
system is a quasi-long range order in phase difference 
when individually the phases are disordered. Such a state 
leads to a counter-intuitive situation of purely topologi- 
cal origin when there are no disspationless individual cur- 
rents while a superflow of oppositely directed two charged 
currents is dissipationless. 

We also stress that the present discussion is relevant 
for much wider class of physical models than two-gap 
superconductivity: the models (1), (2) are also relevant 
in the high energy physics [19, 21]. 
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A note added in proof: We call the reader's at- 
tention to a problem of a BKT transition similar to the 
discussed above type-II limit, which appears in high-Tc 
superconductors and was considered in [29] and refer- 
ences cited therein. We thank G. Blatter for informing 
us about these references. Another, related from a fleld 
theoretic point of view, problem was considered in [30] . 



APPENDIX A 

The pairing of vortices and antivortices of the type 
(1,0) + (0,1) is possible because two condensates are 
coupled by vector potential so gauge invariant equa- 
tions for currents of two types of Cooper pairs are: 

ja = ;^|*a|'V<^a - ^\^a\^A, [a = 1,2) SO that ff 

there are gradients of while there are no gradients of 
4>2 nonetheless coupling by A induces current j2 in the 
condensate *2- For this reason a vortex (1,0) interacts 
with a a vortex with phase winding in different conden- 
sate: (0,1). Let us now consider a vortex (1,0). It has 
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the following particle currents of condensates ^'1,2 around 
its core: 



ji = — i*irv</)i- — 

mi Aruc 



and 
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(28) 



The standard solution for the vector potential of an 
Abrikosov vortex is: 



A=^^|A(r) 



(29) 



where r measures distance from the core and is a unit 
vector pointing along the core. For such a vortex the 
solution for vector potential is 



lAI = 
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For r « A = (e-\/^ + ||) ' eq (30) reads 
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For this reason for the mentioned above vortex pair, in 
a type-II superconductor with the separation r when 
A >> r >> maxj[^j(TBKT)], the interaction mediated 
by the neutral mode due to coupling by vector potential, 
which is composed of oppositely circulating supercurrents 
of two species of Cooper pairs, is depleted because a vor- 
tex (1, 0) at length scale r « A has depleted current j2 
likewise a vortex (0, 1) also has depleted current ji. 
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